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'■^ ' Abstract. We discuss strongly interacting one-dimensional boson-boson mixtures 

C ' and study their quantum phases as the interspecies repulsion is increased. In par- 

O . ticular, we analyze the low-energy quantum emulsion metastable states occurring 

i_~,' at large values of the interspecies interaction, which are expected to prevent the 
system from reaching its true ground state. A significant decrease in the visibility 

rsj , of the atomic clouds is predicted as well, and related to the trapping of the system 

^ ' in the spontaneously disordered quantum emulsion states. 

m 

lO ■ 1 Introduction 

^, 

t^^ ' In recent years the study of systems of ultracold atoms loaded into optical lattices has drawn 

>— ^ , much attention because of their value as experimental probes for the investigation of models of 

condensed matter physics exhibiting key properties of complex many-body systems. One of the 
most relevant examples in this sense is certainly the Bose-Hubbard model, originally introduced 
by Haldane [1] as a simple yet interesting variant of the better known fermionic Hubbard model, 
and later investigated at length and proposed for the description of superfluid ^He trapped in 
porous media in a seminal paper by M. P. A. Fisher et al. [2,. The observation of the distinctive 
C^ ' superfluid-Mott insulator quantum phase transition in a system of ultracold rubidium atoms 

[3] brilliantly demonstrated the prediction that, under suitable conditions, ultracold bosonic 
atoms trapped in optical lattices could provide a physical realization of the Bose-Hubbard 
model [4] . The impressive control in atom cooling and trapping soon allowed the realization of 
more general Bose-Hubbard Hamiltonians, including, e.g., disordered local potentials [5fi\ . 

The quest for novel and unconventional quantum phases has recently pushed the interest 
towards even more general variants of the Hubbard model, involving mixtures of particles 
obeying either the same or different statistics. Beyond their theoretical appeal, these systems 
are relevant to interesting applications such as implementation of disordered systems |7|8| , 
association of dipolar molecules [9] , schemes for quantum computation 10] and realization of 
quantum spin chains and arrays [ll|12j . 

So far, most of the experiments on lattice-atom mixtures have involved atoms obeying 
different statistics. Several experimental realizations of lattice Bose-Fermi mixtures have been 
engineered |8|13|14j . whereas, to the best of our knowledge, the first experiment on mixtures 
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involving two different bosonic species has been carried out very recently at LENS, in Florence 

At the theoretical level, various investigations have been carried out on lattice Bose-Fermi 
[16117118119120121] as well as Bose-Bose mixtures [22I23I24I25I26I27J . This wealth of work demon- 
strates that a comprehensive and systematic study of the phases of atomic lattice mixtures 
requires a significant analytic and synthetic effort. Indeed, one should be aware that even in 
the ideal case of a homogeneous lattice, the Hamiltonian describing a bosonic mixture contains 
five parameters, namely the two intra-species interaction strengths, the inter-species interaction 
strength and one hopping amplitude for each species. Moreover, the possible phases depend on 
many physical parameters such as, e.g., the populations of the two atomic species and their 
commensurability with the lattice size. 

For the above reasons, we confine the present analysis to a rather limited region of the phase 
diagram of a boson-boson mixture, which however encompasses a fair wealth of quantum phases 
effects. Specifically, we consider two different bosonic species loaded in a 1-D homogeneous 
lattice, and vary only the interaction among them. 

We assume that one of the species is soft-core and at unitary filling, whereas the other is 
hard-core and at filling 2/5. Our choice is similar to that considered in Refs. [28|29j . where the 
ground state of a Bose-Fermi mixture on a 1-D lattice is analyzed by means of quantum Monte 
Carlo simulations. Recall indeed that on ID lattices hard-core bosons and spinless fermions 
share many features. While we find a qualitatively similar phase diagram, our analysis differs 
from that in Ref. \2§\ at least in two respects. Firstly, we adopt a less quantive yet more 
computationally affordable approach. 

More interestingly, we focus on the configurations that emerge at strong values of interspecies 
interactions. We find that for sufficiently large interspecies interactions the system supports a 
large number of low energy metastable quantum emulsion states exhibiting glassy features 
despite the complete absence of any source of disorder in the Hamiltonian parameters. Such 
a complex energy landscape for lattice bosonic mixtures was pointed out in Ref. [26] and 
further discussed in Ref. [27] in the presence of the inhomogeneous local potential typical of 
realistic setups. Interesting related results are presented in Ref. [24j based on a multi-orbital 
Gross-Pitaevskii approach, while a phase diagram for metastable states is derived in Ref. [3D] 
for dipolar bosons. We carry the analysis of quantum emulsions further on, discussing the 
arrangement and features of the droplets, as well as the interference pattern obtained after 
free-expansion of the atomic cloud. 

The plan of the paper is as follows: In Sec. [2] we recall the Hamiltonian of the system as well 
as the mean-field approach we adopt. After the introduction of the quantities characterizing 
the quantum phases occurring in the system, our results are presented in Sec. [3J In particular, 
we focus on the behavior of all the aforementioned quantities as functions of the interspecies 
interaction. As we mention, particular attention is devoted to the quantum emulsion states 
characterizing the phase-separated regimes. Interestingly we evidence a phase where the soft- 
core bosons behave as a disordered superfluid featuring a global phase coherence but a spatially 
inhomogeneous density. We conclude with a brief summary of our results. 



2 The System 

The arguments of Ref. [4] can be generalized to a lattice loaded with atoms of two different 
bosonic species. Under suitable conditions, the system is described by the two-flavor Bose- 
Hubbard Hamiltonian 



i / i L 



(1) 
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where the lattice boson operators al j, a/_i, and n^^i = al ^ af^i, create, destroy and count atoms 
of the flavor, or species, / at site i. The parameters Uf and Ui 2 quantify the repulsive interaction 
between atoms of the same or different species, respectively (henceforth intra- and inter- species 
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repulsion). The possibly different hopping amplitude of the two species is quantified by the 
parameters J/. Two-species Bose- Hubbard Hamiltonians similar to the one in Eq. ([T]) have 
been considered previously in several works, possibly referring to different internal states of the 
same bosonic species [23] to spin-1 |31|32|33] or dipolar bosons |9|34j . 

Since we are in a strong interaction regime, we assume that the state of the system \^) has 
a Gutzwiller factorized form. 



m-X{\^.). i'A») = E Jtf(«L)"(«Lri0) 
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where |0) denotes the vacuum state, a/_i|0) = for any i. A time-dependent variational proce- 



dure [25] shows that the complex coefficients Cp^ p2 obey a semiclassical dynamics of the form 
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where the symbols ~ and ~ denote complex conjugation and nearest-neighborhood, respectively, 
and the quantities in Eq. ^ are shorthand notations to simplify Eq. ([3]). 

The ground-state of the system is the lowest-energy normal mode of Eqs. ([3]), which can be 
found as the ground state of the mean-field Hamiltonian 



^ = E ] ^^fAr^f,^ - 1 
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subject to the self consistency constraint a/.j = {'F\afj\'F) = {tpj\afj\ipj). Note indeed that the 
eigenstates of Hamiltonian ([5]) have the form in Eq. ^. A few points are worth observing here. 
First of all, the explicit presence of the (species-specific) chemical potentials /i/, which act as 
Lagrange multipliers fixing the total populations Nf. This is explicitly needed since Hamiltonian 
([5]) does not commute with the total populations Nf = ^ "-/j'l unlike Hamiltonian ([1]). Second, 
the mean-field theory can be equivalently obtained from Eq. ([T]) by means of the so-called 



decoupling approximation, consisting in the substitution al 



/,j«/j 



4,»'^/.i+'^/.'"/j~"/.*"/j- 



More interestingly, we note that the above approach treats the on-site terms of Hamiltonian 
(fT]) exactly. This in particular means that the cross-correlation function Ci^2 introduced in 
Eq. ([7]) below makes sense also for our decoupled trial state ([2]) . We emphasize the explicit site 
dependence of our method. That is, unlike earlier works adopting the Gutzwiller approximation 
[22123136] , we do not force translational invariance on Eq. , and hence do not end up with an 
effective single site theory, which would be unable to describe the structure of phase-separated or 
generically inhomogeneous configurations. Since we are interested in the effects of an increasing 
Ui,2 on a given system, we do not look for the ground state of Hamiltonian ([5]) at fixed values 
of the chemical potentials, as it is usually done. Rather, we adjust the fii and fi2 so that we 
obtain the desired fillings. This procedure, by the way, avoids the so-called species depletion 
problem affecting single-site theories [23136] . 

The quantities ajj, referred to as local order parameters, are directly related to several 
interesting physical quantities, such as the one-body density matrix - related in turn to the 
condensate fraction of the system and to the interference pattern observed in experiments - 
and the superfluid fraction. 
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In the simple case of a single-species homogeneous lattice one recognizes two phases, corre- 
sponding to the local order parameters being all zero or all finite. These situations are readily 
identified with the insulating and superfluid phases, respectively. The resulting phase diagram 
turns out to be qualitatively correct, and its quantitative agreement with the exact result 
improves with increasing dimensionality of the lattice. 

Being an approximation, the mean-field approach suffers from some limitations. For instance, 
site decoupling makes the correlation between two sites of a homogeneous lattice independent of 
their distance, which results in the superfiuid phase being always condensate. This clearly does 
not apply to ID lattices, where the power law decay of the two-site correlations prevents long 
range order in the thermodynamic limit. Despite this artifact, the mean- field approach proves 
useful also on one dimensional systems. Indeed, at the finite — however large — experimentally 
relevant sizes, the mean-field prediction provides an acceptable approximation to the slow decay 
of the exact correlations. This for instance results in qualitatively correct interference patterns, 
characterized by a sharp peak at zero momentum [37] . Such a peak is completely washed out 
in the Mott-insulating regime. 

Furthermore, the local nature of the mean-field order parameters allows for configurations 
where different portions of the system exhibit different phases, thus proving effective in captur- 
ing the geometric fluctuations that add to the quantum ones in the presence of inhomogeneity. 
Different phase domains typically occur in the presence of site-dependent (possibly random) 
local potentials, but are possible also on homogeneous lattices when the translation invariance 
is spontaneously broken. 

The comparison of the results presented here and in Ref. [27] with those in Refs. [26,29] is 
a clear evidence of the validity of the qualitative picture provided by the mean-field approach. 

For the sake of completeness it should be mentioned that the description of a boson-boson 
mixture such that both of the atomic species are in a Mott state reveals a further limitation 
of the mean-field approach. In this case an unrealistic degeneracy for the ground state of the 
system is predicted, which can be resolved only taking into account second-order quantum 
fluctuations [22, 36]. However, the conditions for this artifact to occur are rather specific (the 
total boson population has to be commensurate with the lattice size) and are never considered 
in the following. 



3 Results 

The zero-temperature phases of the system are determined by the properties of the ground state 
of Hamiltonian Eq. |T]). We characterize such properties making use of several quantities. The 
superfiuidity of each species can be estimated as the stiffness under phase variations |20I38I39] 

Fif) = hnr im^lB (6) 

where E{9) is the ground state energy of the Hamiltonian obtained from Eq. ^ by the substi- 
tution at o/.j+i -I- af_ja\ -^-^ — > e'^al -afj^i + e^'^^afjul -_^_^ in the hopping term of species 
/, while Nf — 'Y^Anf^j) is the population of the same species. The introduction of the so-called 
Peierls phases [40' is equivalent to the imposition of twisted boundary conditions. 

Two further quantities will be useful in the characterization of the phases we are going to 
encounter, namely the creation- annihilation cross correlation functions and the local density 
fluctuation 

I,/ i \ i / 

where the operation ~ exchange the two species, i.e. / = 1, 2 if / = 2, 1 respectively. The first 
quantity is specific to mixtures, accounting for nontrivial local quantum correlations between 
particles of different kind and, as we observed in the previous section, is well posed for the 
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Fig. 1. a. Phase diagram of the system as obtained for a lattice with M — 2000. The plotted quantities, 
discussed in Sec.O allow us to identify four different quantum phases; Panels b. and c. refer to regions 
3 and 4 of the phase diagram, respectively, and demonstrate the linear relation between the system 
energy and the number of phase interfaces in the quantum emulsion states the system gets trapped 
into. The panels show 100 data obtained for a lattice with 600 sites. Panel d. shows the color code for 
the phase of the local order parameter characterizing the emulsion droplets in Fig. |2] 



Gutzwiller trial state in Eq. ([2]) • The second quantity is sensitive to spatial inliomogeneity, and 
hence it is useful in the characterization of the quantum emulsion states appearing at large 
interspecies interactions. Also, it conveniently signals phase separation which, in the mean- 
field approximation, occurs through the emergence of ground states breaking the translational 
symmetry of the Hamiltonian due to nonlinear effects. The expected homogeneous ground state 
is recovered as a symmetric superposition of the set of degenerate symmetry-breaking states 
obtained from each other by a lattice translation. 

Finally, the interference pattern in the experimental absorption images of each species is 
basically determined by the Fourier transform of the relevant one-body density matrix summa- 
rizing the first-order correlations [IT] . 
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where Vj is the spatial position of the j'-th optical lattice site, k is the (momentum) coordinate in 
the absorption image. As we will illustrate shortly, these results provide a reasonable description 
of the system, despite the above-mentioned artifact prediction of a finite condensate fraction. 

As we mention above, we focus our attention on a combination of Hamiltonian parameters 
where several different and interesting phases take over as the interspecies repulsion Ui^2 is 
increased. As to the remaining parameters, we choose C/i = 1 as our energy scale, U2 = oo (hard- 
core bosons), Ji = J2 = 0.05, Ni/M — 1 and N2/M = 0.4. In order to rule out contributions 
arising from geometric inhomogeneity we choose a homogeneous lattice comprising M sites, 
with periodic boundary conditions. 

The above described quantities allow us to identify 4 different phases for the system, as it is 
clear from Fig.[TJ When C/1^2 is sufficiently small the situation is similar to the non-interacting 
limit Ui,2 = 0, in which species 2 is supcrfluid and species 1, being strongly interacting and 
at unitary filling, is in a Mott-insulating state. Since the two species are virtually noninteract- 
ing, the cross-correlation Ci^2 vanishes. The different quantities characterizing the system are 
independent of Ui^2 and each species behaves as if it was alone in the lattice. In this situation 
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there are no density fluctuations, as it is expected for single-species systems on homogeneous 
lattices. 

As soon as the interspecies interaction exceeds a first critical value, [/-[ 2 ~ 0.7, it becomes 
sufficiently strong to induce a nonzero cross correlation. The interaction with the hard-core 
bosons melts the Mott insulator, and the atoms of species 1 enter a superfluid phase. This 
can be intuitively explained by observing that, owing to the large interspecies interaction, 
the simultaneous presence of two atoms of different species on the same site is energetically 
unfavourable. Hence the atoms of species 1 are so-to-say expelled from the sites occupied by 
hard-core bosons, and give rise to a nonzero superffuid fraction. Since Ui is still sufficiently larger 
than Ui^2, the bosons of both species remain delocalized and homogeneous, as demonstrated 

by the vanishing of Auf. It is interesting to observe that the increase of Fg in this phase 

(2) 
corresponds to a decrease of Fs , as if the hard-core bosons felt the presence of soft-core 

bosons like a viscous medium slowing down their motion through the lattice. 

When the interspecies interaction exceeds a second critical value, U12 ~ 1.09, the mean-field 
ground-state of the system breaks the translational symmetry of Hamiltonian ([5]), which signals 
the occurrence of phase-separation. This is recognized by the nonzero value of the average local 
density ffuctuations of both species. As we mention above, the expected homogeneous density is 
recovered through a symmctrization of the symmetry-breaking state. As long as the interspecies 
interaction does not exceed the further critical value Ui'2 ~ 1.16, the two atomic species remain 
correlated, C1.2 > 0, and coherent. However, while the soft-core bosons are superffuid, the hard- 
core bosons loose their superfluidity, Fs = 0. For Ui,2 > C/1'2 both Fs and Ci,2 vanish: the 
system has reached a completely phase-separated state where no atoms of species 1 are found 
at sites hosting hard-core bosons, and vice-versa. 

The above scenario is qualitatively very similar to what is obtained in Rcf. 29J , where a Bose- 
Fcrmi mixture is analyzed by means of more precise but much more numerically demanding 
quantum Monte Carlo simulations. The four phases, with a sufficiently large value of U2 instead 
of an infinitely large one, have been proven to be stable under small variations of the interaction 
parameters. 

A very interesting issue recently raised by Roscilde and Cirac [26] concerns the possibil- 
ity that strongly-interacting mixtures exhibit a complex low-energy landscape challenging very 
efficient relaxation dynamics such as that inherent in quantum Monte Carlo simulations. As 
a result, usually efficient minimization algorithms fail to converge to the configuration attain- 
ing the minimum energy, and get virtually trapped into low-energy quantum em,ulsion states 
consisting of a random arrangement of droplets characterized by different phases of the two 
atomic species. The same phenomenology is captured by the Gutzwiller mean-field approach, 
as discussed in Refs. "TT and [30" for strongly interacting bosonic mixtures and single-species 
bosons with long range dipolar interactions, respectively. 

For interspecies interactions exceeding the phase-separation threshold, [/{'j the minimization 
algorithm we adopt gets trapped in quantum emulsion states, which exhibit different features 
in regions 3 and 4 of the phase diagram in Fig.[T] In order to illustrate this, we discuss the prop- 
erties of the quantum emulsion configurations at two representative values of the interspecies 
interaction, Ui,2 = 1.1 and f/1.2 = 1.35. We first of all notice that, as discussed in Refs. [26127] 
the energy of quantum emulsion states exhibits a clear linear dependence on the number of 
interfaces between neighbouring droplets, as it is shown by panels b. and c. of Fig. [T] Clearly, 
the minimum energy is attained by a configuration featuring the lowest possible number of 
phase interfaces, i.e. two. We emphasize that the lowest-energy states in Fig.[T]b., c. have been 
obtained by including an ad hoc constraint on the number of interfaces. Indeed, the presence 
of a large number of low-energy metastable states makes the reaching of the true ground state 
a very hard task. 

Fig. [5] shows the appearance of the typical quantum emulsion state at C/1^2 = 1-1 (left) and 
Ui,2 — 1-35 (right). In both cases the upper and lower panels focus on soft- and hard-core 
bosons, respectively. The solid black line corresponds to the local density {rif.j), whereas the 
coloured areas represent the complex local order parameter a/j = |a/j|e*'^J'-' = (ofj). More 
precisely, the height of such areas is determined by |q;/j p whereas the colors denote the phase 
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Fig. 2. Typical quantum emulsion states in regions 3 (left, f/1.2 = 1.1) and 4 (right, Ui^2 ~ 1.35) 
of the phase diagram in Fig. [T] a. In both cases the top and bottom panels focus on the soft- and 
hard-core bosons, respectively. The black solid line represents the local site occupation (n/j), whereas 
the coloured areas are related to the complex local order parameter, (a/,j) — |a/,j|e"^J'-' . Specifically, 
the height of the area is |a/,jP, while the color is associated to the phase ip/j through the colorbar in 
Fig.[T]d. 



(fifj according to the colorbar in Fig. [T]i. We note that each droplet can be characterized by an 
unique value of the phase, which can vary arbitrarily from droplet to droplet. A light gray solid 
line representing the local density of the species analyzed in detail in the other panel is also 
drawn for comparison in each panel. In both cases one can recognize the presence of "droplets" 
of two kinds. 

For f/1.2 — 1.1 the hard-core species comprises irregularly interleaved superfluid and insulat- 
ing domains, which prevents a superfluid flow. Indeed, on a one dimensional system such as the 

one under investigation, a2,j must be finite at every site in order to produce a nonzero Fs [42]. 
Hence the hard-core bosons are in a kind of glassy insulator, as we discuss in more detail below. 
Conversely, the phase of soft-core bosons could be defined as a disordered superfluid. Indeed 
it also consists of irregularly interleaved domains, which however always exhibit a superfluid 
character since the value of ai^2 never vanishes. 

For C/1,2 — 1.35 the two atomic species never occupy the same lattice site, i.e. they are com- 
pletely phase separated. Mott-like hard-core droplets are randomly intermingled with superfluid 
soft-core domains. Hence in this case the soft-core species is in a Bose-glass-like phase. 

In both the examined cases the droplets composing the quantum emulsion are arranged 
in a random fashion. This feature is quite interesting, and demonstrates the spontaneously 
disordered nature of the low energy configurations in strongly interacting atomic mixtures [26j . 
Recall indeed that the original Hamiltonians do not contain any explicit disorder source. Even 
ignoring the configurations which are equivalent through a lattice translation, a virtually infinite 
number of spatial arrangements of droplets are possible with essentially the same energy, very 
close to the ground state. Furthermore, in cases where a/j is finite only inside one of the two 
classes of droplets — as for species 2 at Ui^2 — 1-1 and for species 1 at Ui^2 — 1-35 — , a further 
source of energy degeneracy lies in the phase (pfj of the order parameter. Note indeed that in 
the above mentioned cases such phase has a constant value within each droplet, but varies across 
different droplets. This is explained by observing that each droplet can be well approximated as 
a virtually isolated and almost uniform system, owing to the vanishing of the order parameter in 
the surrounding droplets. In this situation one expects the chosen periodic boundary conditions 
to be irrelevant. In a uniform system an overall constant phase has no effect on the observables, 
and in particular on the system energy. In the present case different droplets have in general a 
different phase, which manifests itself in observables involving different lattice sites, e.g. like 



^f,jaf,i) 
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(9) 
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This has virtually no effect on the energy, since it gets contributions from terms where t and 
j are nearest neighbours. Indeed when both sites belong to the same droplet the equal phases 
cancel out, while when the sites arc on different sides of a droplet boundary the possible phase 
difference is canceled by the virtual vanishing of the order parameter outside the droplet. It 
is also clear that small variations of the phases tp/^j result in an increase of the energy, which 
demonstrates the local-minimum character of these configurations [43] . 

Interestingly, the different phases of the droplets can play an important role in the exper- 
imental measure of the system coherence. This is usually obtained by imaging the ultracold 
atomic cloud after a few milliseconds of free expansion. The interference pattern in the absorp- 
tion image is substantially described by the Fourier transform of the one-body density matrix 
summarizing the two-site correlations ^4T, Eq. ([5]), whose mean- field form is 



'5'/(k) = ;^Xl("/J~l"/^jl^) 
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M 



^e»(k-r,-..,„)|^^^^.| 



(10) 



We once again recall that in the homogeneous case the mean-field approach results in two- 
site correlations independent of the site distance, which leads to an artifact finite condensate 
fraction in the thermodynamic limit of one dimensional systems. On finite systems this artifact 
is less serious, since the exact result is anyway very slowly decaying with inter-site distance 
[37] . For instance, in both the mean-field and exact cases >5'/(k) exhibits a large peak at k = 
surrounded by small secondary peaks. Therefore the (incipient) long-range order predicted by 
the mean-field approach turns out to be a reasonable approximation to the expected quasi- 
long-range order on homogeneous systems. 

In the case of single a quantum emulsion state such as those depicted in Fig. [2| the correla- 
tions in Eq. ([9|) would be also finite at arbitrary distance, however modulated by a phase factor. 
This is in contrast to what one would expect in a disordered, strongly interacting regime. Nev- 
ertheless, the expected short range correlations are to some extent captured by our mean-field 
approximation. Indeed, it seems reasonable to assume that the state of the system consists of 
a superposition of many almost-degenerate low-energy quantum emulsion states. These will be 
characterized by different phase-interface arrangements compatible with the typical droplet size 
dictated by the system parameters, as well as by different phases of the local order parameter 
within each droplet. This results in pij « |ai||aj| = IctiP > for sites i and j belonging to the 
same droplet, and pij — otherwise, which shows that the correlation lenght is comparable 
with the average droplet size, as expected. When this size is of the order of a few lattice sites 
we expect the absorption image of a quantum emulsion phase to be hardly distinguishable from 
that of a Mott-insulating phase. This presents some analogies with the case of the Bose-glass 
phase induced in a single-species system by (quasi) random local potentials [5'44]. 



4 Conclusion 



In summary, we have discussed the occurrence of four different phases in a ID lattice bosonic 
mixture as a result of the increase of a single parameter, i.e. the interspecies repulsion. The 
features of the different phases have been analyzed by employing suitable observables, and 
special attention has been devoted to the quantum emulsion metastable states occurring at 
large value of the interspecies repulsion [26 27 . Also, we discussed experimental issues such as 
the visibility of these quantum emulsion states, which could greatly reduced by the randomness 
of the phase of the order parameter within each isolated droplet. Owing to the large parameter 
space, we expect that further quantum phases can crop out in different regimes. For example, 
we verified that for a sufficiently large increase in the hopping amplitudes phase 3 is replaced 
by a different quantum emulsion where both species behave as a disordered fluid. Such complex 
phases and the effect of the dimensionality on phase segregation will be the subject of future 
work. 
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